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g -N ' We propose a generic mechanism for the formation of narrow rings in rotating systems. 

For this purpose we use a system of discs rotating about a common center lying well outside 

the discs. A discussion of this system shows that narrow rings occur, if we assume non- 

f^ ' interacting particles. A saddle-center bifurcation is responsible for the relevant appearance 

of elliptic regions in phase space, that will generally assume ring shapes in the synodic frame, 

^\| , which will suffer a precession in the sidereal frame. Finally we discuss possible applications 

^ ■ of this mechanism and find that it may be relevant for planetary rings as well as for semi- 

fT^ ' classical considerations. 

o' 

fvg , §1. Introduction 

o. 

O^ ' With the discovery of ring structures around all mayor planets ^' the phenomenon 

has past from being a breathtaking anomaly in the sky to a rather common feature. 
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While there exists an extended body of work studying features, origin and destiny of 

i-^ '■ planetary rings, ^' ' ^' we believe that there is room for a more general question: How 

generic is the formation of narrow rings? We shall here discuss a model which is quite 

c5 ' remote from any true representation of planetary rings. The purpose will be to show 

i-G ■ that very simple rotating systems can and usually will have a phase space structure 

that allows for the occurrence of rings, even if the particles of such a ring do not 

move on trajectories that could be approximated by the rings shape. Understanding 

the mechanism in a simple model should reveal the genericity we seek. 

We shall start by generalizing the model of rotating discs which we proposed 
earlier. ^' In section 3 we discuss the stability of periodic orbits in this system 
and the ring structures we find as a consequence of stable motion for particles not 
interacting among themselves. Next, we shall see how the introduction of a second 
disc on a smaller orbit may limit the number of rings to as few as a single one. In 
particular we shall see that we do not require both discs to move at the same angular 
velocity to maintain the ring structure. We proceed finally to consider the genericity 
of the results obtained in this particular model. ^^ In particular, we shall discuss the 
consequences of using smooth potentials and even attractive ones such as the 1/r 
potential that occurs in the planetary rings as well as possible applications. 
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§2. The rotating two-disc scattering billiard 

We shall study the planar motion of collisionless point particles in an open 
billiard which consists of two non-overlapping discs on circular orbits (Fig. 1). We 
consider discs of radii di whose center are at distances Ri from the origin (i = 1,2). 
For simplicity, we shall discuss the case where the discs move with the same angular 
velocity cj, which permits a time independent Hamiltonian formulation, although the 
dynamics of the billiard will allow us to relax this condition. The (constant) angle 
formed by the vectors pointing to the centers shall be denoted by (3. In what follows, 
we shall refer to the inner and the outer discs as disc 1 and disc 2, respectively. The 
non-overlapping condition reads as Ri + di < R2 — d2- 




Fig. 1. Geometry of the two rotating discs scattering model: Ri {i — 1,2) define the position of 
the centers of the discs and di their radii; /3 defines the relative positions of the discs. 



In a rotating frame (synodic frame) , the dynamics are given by the Hamiltonian 
1 



J = 2^Px+ Py) - ^{xPy - yPx) + y{x, y) 



2-1 



with V(x,y)=0 outside the discs, and infinite inside. The explicit time independence 
of J in Eq. 2-1 implies that it is a conserved quantity of the global flow, which 
is known as the Jacobi integral. We shall define the coordinates such that in this 
rotating frame the outer disc lies on the positive x-axis. We note that, if the discs 
move with different angular velocities, there is no first integral of motion. 

On a fixed frame (sidereal frame), the dynamics are straightforward. A particle 
moves on a straight line with constant velocity until it hits one disc, which changes 
both the magnitude and the direction of the velocity vector. ^' After the collision, 
a point particle wins (looses) energy if it bounces off the discs on their front (back) 
side. This naive observation is precisely the key point for the construction of the 
symmetric periodic orbits of the system. If a particle bounces radially, that is at one 
of the intersection points of the disc and the line that joins the center of rotation 
and the center of the disc, then the energy is conserved and incoming collision angles 
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are equal to outgoing ones. If we have initial conditions which display consecutive 
radial collisions with the disc, they belong to a trapped orbit. In the synodic frame, 
such trapped orbits are actually periodic orbits which preserve the symmetry of the 
problem. ^' As shown in Ref. 8), the fundamental periodic orbits in the construction 
of the chaotic saddle are precisely the symmetric periodic orbits, or deformations of 
them for the non-symmetric cases. 

We shall first consider either of the two discs by itself. In terms of the Jacobi 
integral Ja [i = 1,2), their symmetric periodic orbits are given by^-* 



^h 



uj'^{Ri - d, 



\2 



[{2n + 1)7T - 2aif 



2cos^ai-{{2n+l)TT-2ai)sm2ai . (2-2) 



Here, Oj G [—it/2, tt/2] is the outgoing angle for the radial collision with the disc 
(negative values correspond to a particle running against the rotation), and n = 
0, 1, 2, . . . denotes the number of full turns that the disc completes before the next 
collision. These curves are plotted for some values of n in Fig. 2a as continuous 
lines. Notice that n defines continuous families of such periodic orbits with a precise 
hierarchical arrangement. 

Similar families of periodic orbits exist which alternate bounces among the discs. 
However, symmetry considerations of these two-bounce periodic orbits imply that 
they retain the symmetric character, i.e. bounce radially, only for /3 = and /3 = vr. 
That is, except for those values of /3, related periodic orbits are alternatively accel- 
erated and decelerated after collisions with the discs. The two-bounce symmetric 
periodic orbits are given by 

Ji2 = -Vp - oj{R2 - d2)vp sin 02 (2-3) 

where the (constant) velocity between collisions is (for /? = 0, vr) 

a;(i?2 — ^2) sin(ai — 02) 



[niT + j3 + ai — 02] sin ai 

In these cases ai and 02 are geometrically related by the identity 

sinQ2 sincti 

Ri+di~ R2-d2' 



(2-4) 



(2-5) 



Notice that Eq. (2-5) defines the limits \a2\ < a™^^ for which the orbits J12 
exist, where sin(a™'^^) = (i?i + (ii)/(i?2 — c?2)- The case where the equality holds 
correspond to tangent collisions with disc 1. Again, n is defined as the number of 
turns completed by the discs during one period of the symmetric periodic orbit. 

In Fig. 2b we have plotted the symmetric periodic orbits which involve disc 2 
for /3 = for some values of n, as given directly by Eqs. 2-2 and (2-3-2-4). It follows 
from Eqs. (2-3-2-4) that for (3 = tt the n curves correspond to the n + 1 curves of 
(3 = 0. As it can be appreciated in Fig. 2b, the curves J12 display cusps at the edges 
of the interval where they are defined (|a2| = af"^^) and decrease monotonically. 

We notice that the curves J22 (continuous lines in Fig. 2a) are convex or concave 
depending only on the sign of 02, and then display local maxima and minima. The 
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Fig. 2. (a) Periodic orbits J22 (continuous curves) and J12 for /3 = (dotted curves) for various 
values of n (-Ri = 1.8, di = 0.1, -R2 = 3, ^2 = 1)- Notice tiie divergence of J12 as a — > for 
the n = case, (b) Enlargement of (a) showing the periodic orbits of J12 for j3 = (continuous 
curves) for various values of n; some periodic orbits J22 are shown as dotted curves. The family 
of periodic orbits for l3 — n are obtained from those for /3 = by rearranging the labels (see 
text) . 



important fact to observe here is that the derivative of the J12 famihes is strictly 
negative for any n, while dJii/da2 (i = 1,2) vanishes at least in one point for all 
n. This vanishing derivative is important in the context of stability of the periodic 
orbits and their bifurcations. 



Narrow rings 5 

§3. Stability of the J22 symmetric periodic orbits: formation of rings 

In order to consider the stabiUty of the periodic orbits, we shall begin by ana- 
lyzing those of the outer disc. As noticed by Henon^-* for the restricted three-body 
problem, the stability properties are related to the value of dJu/dai. In fact, a 
saddle-center bifurcation scenario associated with the maxima and minima of the 
curves Ju is responsible for the appearance or destruction of the symmetric peri- 
odic orbits for the rotating one-disc billiard. ^" ^^' Therefore, just below each of the 
maxima or above each of the minima, there exists an interval of Jacobi integrals 
where one of the corresponding periodic orbits is stable and surrounded by invariant 
tori; the other fixed point corresponds to a hyperbolic one whose manifolds build 
an incomplete horseshoe. These maxima and minima are defined by the condition 
dJ22/da2 = 0, which can be rewritten as 



tana 



(2n + l)7r - 2a* 



±1. 



(3-1) 



This equation defines for a given n the angles 02 = a* where the maxima or minima 
are found, independently of the specific geometrical parameters considered in the 
model. 




0X102 0,00: 



Fig. 3. Poincare surface of section showing a stability island (J = 1.173). An ensemble of initial 
conditions started within this stable region will not escape along scattering trajectories. 



In Fig. 3 we present the surface of section for the one-disc rotating scattering 
billiard for J22 = 1.173, illustrating the regions of stable motion just described. We 
have chosen the parameters R2 = 3 and 1^2 = 1 (lj = 1). The surface of section is 
defined here by the angles a and (j). The latter angle characterizes the position of 
the collision point on the disc. ^' There are only two fixed points of period one for 
this value of the Jacobi integral, which correspond to the periodic orbits near the 
maximum of the n = case (see Fig. 2a). One of these fixed points is stable and 
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appears surrounded by KAM tori. These KAM tori as well as some regions of chaotic 
motion are bounded dynamically by the invariant manifolds of the partner hyperbolic 
fixed point. This implies that, even-though the problem is actually a scattering 
billiard system with a (hard disc) repulsive potential, there exist initial conditions 
of stable motion which are trapped by successive collisions with the rotating disc. 

The KAM tori shown in Fig. 3 are robust and persist under small changes of 
the Jacobi integral. For larger changes, bifurcations occur but those of the central 
elliptic point are at first hidden from scattering trajectories by surviving KAM sur- 
faces. This scenario persists under variations of the Jacobi integral until the formal 
parameter describing the (local) binary horseshoes development attains the value 
2/3 (see Refs. 11), 12) for details), which indeed defines a continuous interval for 
J where bounded trajectories exist. As we further change the value of the Jacobi 
integral we reach stages of the incomplete horseshoe that correspond to hyperbolic 
structures and others where very small stable islands reappear. 

In the present context, we shall concentrate on a set of initial conditions that 
remain trapped under small perturbations. That is, we consider all those initial 
conditions which belong to a bounded region, without considering if their actual 
motion is periodic, quasi-periodic or chaotic. This set of initial condition has strictly 
positive measure, i.e. there is a non-zero probability for generic initial conditions 
of being within this set-^^^. Note that this situation is essentially different from the 
hyperbolic one, where at most a Cantor set of measure zero is trapped whether the 
horseshoe is complete or incompleted^'. 

We shall stress that, for a generic ensemble of initial conditions defined on the 
whole phase space, which is not restricted to a fixed Jacobi integral shell, most ini- 
tial conditions will escape to infinity. However, particles whose initial conditions are 
chosen close enough to the stable periodic orbits and within the bounded regions 
described above, will be trapped by collisions with the disc. This observation is fun- 
damental for understanding the occurrence of patterns or rings formed by infinitely 
many particles whose initial conditions are started precisely on dynamically bounded 
regions found in distinct intervals of the Jacobi integral. 

In Fig. 4a we present in the sidereal frame the ring formed by an ensemble of 
initial conditions whose Jacobi integral belongs to a neighborhood of the minima 
of J22 for n = 1. Figure 4b displays the set of rings formed by the minima of J22 
for n = 1,2. In these figures, each point represents at a given time the position 
of a single particle which moves on a rectilinear trajectory. These particles have 
collisions with the disc at different times and in any point along the disc's circular 
orbit (sidereal frame). The loop shown in Figs. 4b is associated with the n = 2 ring. 
It results from particles on different rectilinear trajectories (with approximately the 
same Jacobi integral) which come closer during some interval of time. In general, 
there will be n loops associated with the retrograde (02 < 0) J22 symmetric periodic 
orbits, while there will be n— 1 loops for the corresponding prograde orbits (02 > 0). 
The isolated points in both figures correspond to initial conditions, that were not 
located in the bounded regions, but have not escaped after the time at which the 
computation was cut. Another remarkable aspect of Figs. 4 is that instead of viewing 
it as a stroboscopic map of many non-interacting particles in the sidereal frame, we 
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Fig. 4. Examples of rings in the sidereal frame associated with the outer disc (shaded regions) for 
the rotating two-disc scattering problem. The parameters of the inner disc are defined to filter 
all but the rings shown (the parameters of the outer disc are the same as in Fig. 2). (a) Ring 
corresponding to the n — 1 minimum of the J22 curves {Ri = 1, di = 0.52); (b) corresponding 
rings for the n = 1,2 minima (_Ri — 0.96, di — 0.52). The dashed circles correspond to the 
inner and outer circles defined by the orbits of the outer and inner discs, respectively. 



could think of the narrow bands of points as curves, and would obtain a picture 
of the elliptic central orbits of the stable islands of a single particle in the synodic 
frame. 

We can ask about possible structure in the rings of our model. First of all, 
different rings may coexist, as we see in Fig. 4b. Their particles move at different 
speeds, and thus probably are not good candidates for the braids of planetary rings. 
However, we could expect them to be important, say in a semi-classical context, 
where we might find interference between two structures. Another type of structure is 
related to the period-doubling scenario, that happens as the elliptic orbit bifurcates. 
There, several strands of particles may be intertwined, while they move roughly at 
the same speed and maybe confined to a single ring by surviving exterior KAM 
surfaces. The details of the structures found may be quite complicated and will 
depend on the dynamics. Yet the very fact that such structures generically appear 
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is significant. 

§4. Stability of the rings for the two-disc model 

It is important to stress that the appearance of rings as described above re- 
quires one rotating disc only. Indeed infinitely many rings exist with ever increasing 
numbers of loops, which give rise to arbitrarily complicated structures in wide rings. 
For the full two disc rotating billiard such rings may appear associated with either 
one of the discs. In what follows we shall not consider the rings associated with 
the inner disc only, as we are interested in narrow rings between the two discs that 
are actually influenced by both of them. We shall thus concentrate on the periodic 
orbits associated with the outer disc, which may or may not have collisions with the 
inner one. The characteristic curves for the J22 and J12 symmetric periodic orbits 
are shown in Fig. 5 for angles /3 = 0, vr between the discs; the size effects of the inner 
disc are taken into account (i?i = 1.8, di = 0.1). 

As mentioned in Section 2, the J12 (Fig. 2a) characteristic curves for families 
that bounce on both discs decrease monotonically as a function of 02. In fact, it can 
be shown that dJi2/da2 is strictly non-zero. Therefore, these families of periodic 
orbits do not undergo any bifurcation as a function of the Jacobi integral. Yet it 
was near the bifurcations that we found stable islands for the one-disc orbits. Indeed 
we do not find any for the two-disc orbits J12 and, as we shall see in the next 
section, there is no reason to expect any. By consequence, these families of periodic 
orbits do not contribute to the formation of rings. Nevertheless, their presence is 
intimately related to the disappearance of orbits of the J22 type, as can readily 
be inferred from Fig. 5. This erosion due to the inner disc will in general eliminate 
some of the extremal points of the characteristic curves; if a sufficiently large interval 
near the extrema is cut the stable island will no longer exist. To see this, consider 
initial conditions belonging to some symmetric periodic orbits J22. Depending on 
the geometrical parameters of inner disc, some of these orbits may collide with the 
inner disc. The cusps displayed by the J12 families indicate the tangent collisions 
with the inner disc (see Figs. 2 and 5). The particles that collide with the inner disc 
will loose the needed correlations to build a symmetric periodic orbit, so after a few 
more collisions these particles will typically escape. 

A gap in the characteristic curves for the J22 orbits thus occurs. Its location will 
depend on the angle and radial distance at which the two discs are positioned, while 
its size will depend on the radius of the inner disc. To illustrate this more clearly 
Fig. 6 shows the J22 curves only, for some angles (3. 

If the inner disc moves at a speed that is incommensurable with the outer one 
we find a vastly simplified picture: It will sweep all periodic orbits that penetrate the 
circle described by the outer edge of this disc. This erosion mechanism is governed 
in our model by q™^^ compared to the solutions of Eq. 3T (which are independent of 
the geometrical parameters of the model). Its typical time scale depends on the size 
of the inner disc and the time dependence of /?. The absolute value of the solutions 
of Eq. (3-1) converges to 7r/4 as n ^ 00, from above for prograde and from below 
for the retrograde orbits. This arrangement of the solutions implies that the lowest 
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Fig. 5. Symmetric periodic orbits for the two-disc model including the finite size effects of the inner 
disc, (a) /3 = 0; (b) /3 = tt. 



prograde solutions will be the first to survive, if we start with very small radial 
distance of the discs and gradually increase it. ^-* Note though that the n = curve 
in Figs. 2 has no stable prograde periodic orbits in this model (it corresponds to a 
whispering gallery orbit). If we want to restrict our considerations to a few rings we 
are forced to look at the prograde part of the branches for the lowest n starting from 
n = 1. Indeed, in Fig. 4a we see a single prograde ring with a circle drawn exactly 
in the path of the outer edge of the inner disc that ejects all particles that might be 
associated with other values of n or with retrograde motion. In Fig. 4b the prograde 
rings for n = 1 and n = 2 survive. These figures where generated by the dynamics 
with two discs present and somewhat scattered initial conditions. The points that 
are not on the rings had unstable initial conditions and are on outgoing scattering 
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Fig. 6. Periodic orbits of the J22 family (n = 0, 1...4) for different values of /?. (a) /3 = 5°, (b) /3 = 
10°, (c) 13 = 30°, (d) /? = 45°, (e) /3 = 90° and (f) /3 = 135°. 



orbits, but have not yet left the domain of the figure for the time we integrated the 
equations of motion. 

Note that for different frequencies of the two discs the Jacobi integral is no longer 
an invariant of the entire system, but in the subspace of phase space where our rings 
live the integral is still conserved. 
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§5. Genericity of the scheme and its importance for planetary rings and 

other appHcations 

If we consider Fig. 2a, we see that symmetric periodic orbits of the J22 family 
only exist between an upper and a lower limit for the Jacobi integral J. The stable 
islands discussed in the previous sections are a consequence of the generic scenario 
of saddle-center bifurcations. Thus we expect and indeed find at each branch point 
an elliptic an a hyperbolic periodic orbit. The elliptic one will then turn into an 
inverse hyperbolic fixed point via the period doubling route. 

How generic is this scenario? The presence of the cusp where two-disc orbits and 
one disc orbits touch in their characteristic curves is certainly non-generic. Thus the 
instability of the two-disc orbits and the behaviour near that cusp will have to be 
discussed in each particular case with care. As long as the cusp is not too near the 
maximum of the generic curves, the saddle-center bifurcation will proceed generically 
Cases where the two are close will need special attention. 

If we consider potential mountains instead of hard discs, we will find an alter- 
native generic scenario; the sudden bifurcation ^^' associated with the top of the 
mountains. Yet, the scenario we discuss will still exist for trajectories that do never 
reach the region of phase space where the sudden bifurcation occurs. 

The question becomes more involved if we consider attractive potentials. In view 
of the origin of this problem in planetary rings we shall focus on the 1/r case. We 
know from Henon's work ^^^ and some of our own ^^' that the so called consecutive 
collision orbits govern the entire chaotic saddle for the restricted three-body problem 
with small mass parameter /x, i.e. small ratio between the masses of the lighter and 
heavier binary components. These collision orbits are a good approximation both 
for the 1/r problem with vanishing mass parameter and for the hard disc problem 
with vanishing disc radius. In fact, they coincide for vanishing disc radius, when we 
consider an overall central 1/r attractive potential. Therefore, the two problems are 
not quite as distinct as one might believe. 

The main diff'erence is that these so called collision orbits may extend both 
outward and inward from a moonlet. If we consider the orbits that extend inward 
from the outer moonlet, we are in a situation quite analogous to the one we found 
for the discs, and the inner moonlet will sweep the ones that reach inside its orbit. 
If this last observation is to be effective in finite time, we have to pass to finite but 
small /i. For the orbits that reach outward from a moonlet, the two will interchange 
their role. The outer one will sweep the collision orbits of the inner moonlet that 
reach outward beyond the orbit of the outer one. In Fig. 7 we illustrate, in the 
sidereal frame, both situations with one collision orbit that intersects the orbit of 
the other moonlet and another that does not. Clearly the zero mass approximation 
gives no information about stability. Calculations to shed light on this aspect are 
under way. 

Thus, despite of the caveats we pronounced initially, the mechanism we discov- 
ered may well be of relevance for the structure of planetary rings. In particular due 
to the generic nature of the elliptic regions, these are expected to be more stable 
against perturbations than the orbits of a super integrable system, and could thus be 
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Fig. 7. Examples of consecutive collision orbits in the sidereal frame generated by collisions with 
(a) the outer moonlet and, (b) the inner one. The dotted and dashed circles represent the 
circular trajectories of the moonlets. Note that in each picture one of the orbits will be swept 
by the other moonlet. 



associated to a prediction that older ring structures would tend to consist of narrow 
rings. 

Returning to the more general aspects, the genericity of the occurrence, which 
we suspect also for attractive potentials, could imply importance in other systems 
such as rotating molecules or high spin states in deformed nuclei in the framework of 
semi-classical theory. In this case, the mean field potential can be assumed to rotate 
stiffly and the survival of orbits may depend on the detailed geometry we face. To 
simplify matters in this situation, the Jacobi integral is a constant of the motion. 
This may be of interest when the adiabatic or Born-Oppenheimer approximation 
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breaks down. In nuclei this can well occur for high spin states; the deviations of 
level statistics from the compound nucleus prediction observed for such systems ^^' 
was attributed to a harmonic component simulated by the centrifugal force, but it 
is worthwhile to consider the effect of integrable islands of the type mentioned. In 
molecules it seems less obvious to find a place for non-adiabatic effects, yet Rydberg 
molecules or molecules with exotic leptons may provide interesting examples. 
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